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ABSTRACT 

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of 
currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, 
hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively 
well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence 
of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak 
magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus 
peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters 
that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the 
brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The 
dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no 
evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P. 

Key words. Stars: supernovae: individual: LSQ13cuw 


1. Introduction 

spite of the surge in supernova (SN) discoveries over recent 
ars, a number of persistent issues have remained unsolved for 
decades. In order to gain insights from even larger samples of SN 
discoveries that are expected with the next generation of tran¬ 
sient surveys, alternative approaches may be necessary. Despite 
significant expenditure of effort, one issue remains: the lack of a 
complete and consistent mapping between the evolutionary stage 
during which a massive star explodes, and the type of SN that is 
observed as a result. 

Most - if not all - SNe can be classified as belonging to 
the Type II or Type I categories depending, respectively, on the 
presence or absence of hydrogen in their spectra. A further set 
of four (“a”,“b”, “c”, “n”) spectroscopically motivated classifi- 
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cations provides a grouping into distinct subtypes to denote the 
presence or absence of either hydrogen or helium, or both. The 
“n” indicates the presence of narrow (< lOOOkms -1 ) hydrogen 
lines in emission, and is equally applicable to Type I and II SNe. 
Although assigning a particular classification to rare peculiar 
SNe or transitional objects may not be trivial, the vast major¬ 
ity of SNe can be typed on the basis of a single spectrum. It has 
long been recognized ( iBarbon et all 1979!) that two further divi¬ 
sions of Type II SNe can be made, based solely on light curve 
morphology post maximum brightness; these are the II-plateau 
(P) and Il-linear (L) varieties. As their names imply, the former 
displays an extended period (typically ~100d) during which the 
(VRZ-band) brightness remains constant, usually to within < 0.5 
magnitudes; Type II-L SNe decline linearly at a rate of a few 
hundredths of a magnitude per day, before settling onto the ra¬ 
dioactive tail powered phase. 
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In the canonical picture (e.g. iNomoto et al.lll995h . the ob¬ 
served diversity in the spectroscopic and photometric behaviour 
of core-collapse SNe is explained as a one-parameter sequence 
(II-P —> II-L —> lib —> lb —> Ic) governed primarily by the hy¬ 
drogen envelope mass of the progenitor at the time of explo¬ 
sion. Thus, in a single-star scenario, the expectation would be 
for higher mass stars to evolve to a stage with the thinnest of 
hydrogen envelopes. This picture is intuitively appealing, and is 
borne out by observations of the progenitors of several Type II-P 
SNe, which have been shown to have red supergiant progenitors, 
with masses in the range of 8.5 < M/M Q < 16 (ISmartt et al.1 
120091 and references therein). Barring some exceptions, direct 
detections of the progenitors of other subtypes of core-collapse 
SNe have remained elusive. 

For stars with massive hydrogen envelopes, the plateau phase 
is attributed to the balance between cooling and recombination in 
the ejecta dGrassberg et al.l[T97 ll) . A recombination wave sets in 
when the ejecta have cooled to -5-6 kK, the recombination tem¬ 
perature for hydrogen. The cooling/recombination wave moves 
in a manner such that a roughly constant amount of internal en¬ 
ergy is advected through it per unit time, resulting in a constant 
luminosity. When the entire envelope has been traversed, the 
above no longer holds, marking the end of the pleateau phase. 

Suppressing the cooling/rec omb i natio n w_ave results_in a 
linearly declining light curve dBlinnikov & Bartunovi 1 19931) 
and can be achieved if ejecta densities are sufficiently low 
dGrassberg et al.llT971 ). A natural way of fulfilling this require¬ 
ment is to reduce the hydrogen envelope mass to no more 
than a few solar masses as has been shown by seve ral au¬ 
thors (S wartz et al. 119911 Blinnikov & Bartunovm993 ). How¬ 
ever, large progenitor radii are necessary in order to attain the 
brightness required by observed light curves. 

Thus, within the framework summarized above, one would 
expect Type II-L SNe to arise from stars that have lost more of 
their hydrogen envelope than the progenitors of Type II-P SNe, 
thereby giving rise to linearly declining light curves. Further¬ 
more, one might expect there to be a continuum of decline rates 
between these two groups of Type II SNe. 

Over the decades, several_studies have considered this very 
issue. Earlier work bv lPatat et al.l (1 19941) . who considered a het- 
eregeneous sample of Type II SNe, found the /I-hand light curves 
of Type II-P and II-L SNe at epochs of <100d to overlap in 
terms of their decline rates. Since then, larger datasets of well- 
sampled lightcurves have become available. However, complete 
consensus on whether the Type II-P and II-L SNe show a truly 
bimodal distribution in decline rates is not e vide nt. The larger 
(Type II) samples of I Anderson et al.l d2014bl) and I Sanders et alJ 
( 20141) seem to in dica te a continuum of decline rates, while 
lArcavi et al.l (120121) and lFaran et al.l (12014bh suggest that distinct 
subtypes may be apparent if a judicious choice of decline rate 
post-peak brightness is made (> 0.3 - 0.5 mag , decline at -50 d) 
i.e., similar to that proposed bv lLi et akl (i201 ll) . 

The evolution of the light curve from explosion to peak 
brightness contains information about the size, and therefore the 
type of star that exploded. During the early radiation-dominated 
phase, the luminosity is d i rectly proportional to the pro genitor 
radius ( Swartz et al.lll99lt iBlinnikov & Bartunovlll993l) . Thus, 
considering the shape and the total time of the rise to peak bright¬ 
ness for Type II-P and -L SNe, is a potentially fruitful endeavour. 
Indeed, for thermonuclear SNe, the exceptionally early discov¬ 
ery of SN 201 lfe, led to the firs t ob ser vatio nal confirmation of a 
white-dwarf progenitor dNugent et al.ll201 lb . 

For almost all types of SNe, such observations are particu¬ 
larly challenging given the combination of faintness at the time 


of explosion, and the often rapid rise to peak brightness. Still, the 
type II-P/L SNe are most amenable to such studies, as the 56 Ni- 
powered phase sets in later than in type I SNe. Nevertheless, the 
rise time behaviour of Type II-L SNe is largely uncharted terri¬ 
tory. 

Motivated by the discovery of LSQ13cuw soon after explo¬ 
sion, we embarked upon a study of Type II-P/L SNe with well- 
constrained explosion epochs. In what follows, we consider the 
photometric and spectroscopic evolution of LSQ13cuw together 
with a carefully selected sample of SNe assembled from the lit¬ 
erature. In order to identify global trends, we compare the obser¬ 
vations with two well-known and physically intuitive analytical 
models. 

A study addressing si milar issues was recently presented 
bv lGonzalez-Gaitan et al.l (120151) . Although we refrain from de¬ 
tailed comparisons which are beyond the scope of this paper, we 
note that in spite of broad agreement, there are key differences 
in the selection criteria of the SN sample, the approach, and the 
interpretation. We endeavour to highlight the most significant of 
these in what follows. 

The paper is divided into three relatively self-contained 
parts: LSQ13cuw observations are presented in ^2 the full sam¬ 
ple of SNe and comparison with models forms while a fuller 
discussion of the analytical light curve models and the influence 
of key parameters appear in the Appendix. 


2. Observations and data reduction 

Supernova LSQ13cuw was discovered on 2013 October 30, dur¬ 
ing the course of the La Silla Quest Supernova Survey (LSQ 
iBaltav et al J 120131) . The SN was undetected down to a limit¬ 
ing /-band magnitude of 21.21, two days prior to the discov- 
ery. It was spectros copically classified by the PES STO survey!! 
( Smartt et al.ll2014l ) first as a Type I SN at z -0.25 dBersier et ak 
1201 3l) given the combination of a featureless blue continuum, 
and a long rise time. About a week later, based on second spec¬ 
trum, it w as r e class ified as a Type II-P SN at maximum at 
z -0.05 ( iPennefeld et al.l2013l) . Continued monitoring revealed 
a type II-L light curve ( 32.31) . As part of the PESSTO follow¬ 
up campaign, we obtained 5 spectra ranging from 25 to 84 days 
after explosion; additional imaging up to -100 days after explo¬ 
sion was acquired from a number of different facilities. 

2.1. Data reduction 

Optical photometry was obtained primarily with the Optical 
Wide Field Camera, 10:0, mounted on the 2m Liverpool Tele¬ 
scope (LT; g'r'i' filters), with additional epochs from the Las 
Cumbres Observatory Global Telescope Network lm telescope 
(LCOGT; g'r'i' filters) and the ESO Faint Object Spectro¬ 
graph and Camera, EFOSC2, mounted on the 3.58m New 
Technology Telescope (NTT; V#641, g# 782, r#784, (#705 fil¬ 
ters). All data were reduced in the standard fashion using the 
LT/PESSTC0/LCOGT pipelines, including trimming, bias sub¬ 
traction, and flat-fielding. When nece ssary, cosmic rays were re¬ 
moved using the lacosmic algorithm (Ivan Dokkumll200ll) . 

Point-spread function (PSF) fitting photometry of 
LSQ13cuw was carried out on all images using the cus- 


1 Public ESO Spectroscopic Survey of Transient Objects; 
www.pessto.org 

2 The PESSTO pipeline, developed by S. Valenti, comprises a set of 
python scripts which call pyraf tasks to reduce EFOSC2 data. 
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tom built SNOoP\Q package within iraiQ. For the g'r'i' -bands, 
a number of local sequence stars (see Table lB.ll in the Appendix) 
from the SDSS DR9 cataloging were used to determine colour 
terms on photometric nights and zeropoints in order to calibrate 
the photometry of the SN to the SDSS system. For the NTT 
V-band, stars within standard Landolt fields (lLandolti 1 1992b 
were used for calibration. We estimated the uncertainties of 
the PSF-fitting via artificial star experiments. An artificial 
star of the same magnitude as the SN was placed close to the 
position of the SN. The magnitude was measured, and the 
process was repeated for several positions around the SN. The 
standard deviation of the magnitudes of the artificial star were 
combined in quadrature with the uncertainty of the PSF-fit and 
the uncertainty of the photometric zeropoint to give the final 
uncertainty of the magnitude of the SN. We note that the NTT 
filters are not standard Sloan filters. V#641 and r#784 are Bessel 
filters, while g #782 and i #705 are Gunn filters. These usually 
differ by about 0.1 to 0.2 mag from the corresponding Sloan 
filters. We have not attempted a correction as the difference is of 
the same order of magnitude as the error on the measurement, 
and affects only one epoch in our light curve. 

In addition we analysed pre- and post-explosion images for 
LSQ, which were taken by the ESO 1.0m Schmidt Telescope. 
The same PSF-fitting technique as described above was per¬ 
formed on th e imag es. The images were taken with a wideband 
filter dBaltav et al.ll2013h . and we calibrated the magnitude of 
the SN to SDSS r' by determining zeropoints using a local se¬ 
quence of stars (Table IbTTT ) . However, it was not possible to ap¬ 
ply a colour transformation since there are no available images 
of the SN with other filters. Despite this, it is likely that the rela¬ 
tively large uncertainties of the LSQ magnitudes are greater than 
any differences a colour transformation would make. In the pre¬ 
explosion images we determined limiting magnitudes by adding 
artificial sources to the images at the position of the SN, which 
were created by building a PSF from nearby stars. The mag¬ 
nitude of the PSF was measured and decreased until it was no 
longer detectable (to 3 cr). The magnitude of the faintest (3 cr) 
detection was taken to be the limiting magnitude of the image. 
We found the limits to be +21.21 and +21.19 mag for images 
taken -0.11 d and -2.02 d before explosion. 

LSQ13cuw was independently obs erved by the Catalina 
Real-time Transient Survey (CRTS, iDrake et aDl2009 1 with the 
Catalina Sky Survey (CSS) Schmidt telescope under the Tran¬ 
sient ID CSS131110:023957-083124. We have included the re¬ 
ported unfiltered magnitudes in Table Q] which summarizes the 
photometric observations for LSQ13cuw. The CRTS data are in 
good agreement with our LSQ photometry. 

A series of five optical spectra were obtained with the 
NTT+EFOSC2 (see Table 0. The spectra were reduced with 
the PESSTO pipeline using standard techniques. These included 
trimming, bias subtraction, flat-fielding, wavelength calibration 
via arc lamps, and flux calibration via spectrophotometric stan¬ 
dard stars. The spectra were additionally corrected for telluric 
absorption by using a model spectrum of the telluric bands. 

Features attributable to the Na i doublet from interstellar gas 
either in the Milky Way, or the host galaxy, are not apparent 
in the spectra of LSQ13cuw (see Figure 0. In order to de- 

3 SuperNOva PhotometrY, a package for SN photometry implemented 
in IRAF by E. Cappellaro; http://sngroup.oapd.inaf.it/snoopy.html 

4 Image Reduction and Analysis Facility, distributed by the National 
Optical Astronomy Observatories, which are operated by the Associa¬ 
tion of Universities for Research in Astronomy, Inc, under contract to 
the National Science Foundation. 

5 www.sdss3.org 


rive an upper limit to the equivalent width (EW) of a putative 
Na i D absorption feature, we constructed a weighted stack of all 
LSQ13cuw spectra, where the weights reflect the signal-to-noise 
ratio in the 5400-5700 A region. We then created a series of ar¬ 
tificial Gaussian profiles centred at 5893 A and subtracted these 
from the stacked spectrum. The FWHM of the Gaussian pro¬ 
files was set to 18 A, which corresponds to our lowest resolution 
spectrum. We increased the EW of the artificial profile in steps 
of 0.1 A in order to determine the limiting EW at which an artifi¬ 
cial line profile as described above would be just detectable in the 
stacked spectrum. Using this method we derive the upper limit 
for the EW of the NaiD 4/15890,5896 blend to be < 0.9. Us¬ 
ing the relation in Equation 9 of lPoznanski et al.l (12012h . we find 
E(B - V)host < 0.16 mag. Applying the relation of Turatto et al.l 
(120031) we fin d E(B - V)h ns t . < 0-13 mag. 

Recently, [Phillips et al.l (1201 3l) argued that the EW of the 
NaiD absorption is unreliable as a measure of the dust extinc¬ 
tion from the host galaxies of Type la SNe. Nevertheless, a weak 
or undetectable NaiD absorption seems to be in agreement with 
little or no extinction. 

Other circumstantial factors point towards negligible host 
galaxy reddening. LSQ13cuw lies in an apparently isolated re¬ 
gion at the edge of its presumed host galaxy which is extremely 
faint (see Section 12,2b . It is unlikely that this region is a pocket 
of high extinction in an otherwise unremarkable region. Fur - 
thermore, recent studies, albeit on SNe la (e.g. iHolwerda et al.l 
12015 1. have shown that the Ay is positively correlated with the 
radial position of the SN within the host galaxy. With all possible 
factors taken into account, we assume no host galaxy extinction 
for LSQ13cuw, and ado pt a value of £_(fi - V) = 0.023 mag for 
the Galactic extinction lSchlaflv & Finkbeinej (1201 111 . 

We obtained a redshift of 0.0453 ± 0.0028 for LSQ13cuw 
from the Balmer lines present in our highest S/N spectrum 
(+32d). Assuming Hq = 72kms~ 1 Mpc~ 1 this results in a dis¬ 
tance of 189 Mpc ± 12Mpc and an absolute r'-band peak mag¬ 
nitude of -18.04 ± 0.17. 


2.2. Host galaxy 

LSQ13cuw exploded 13703 W and 0762 S of SDSS J023958.22- 
083123.5 (Figure Q}. In order to determine whether SDSS 
J023958.22-083123.5 was indeed the host galaxy of LSQ13cuw, 
we included it in the slit with the SN on the final epoch (2014 
Jan. 24) of our spectroscopic observations. The spectrum was 
reduced in the same way as the SN spectra (see Figure 0 for 
the reduced and calibrated spectrum). It is almost featureless, 
so a reliable redshift cannot be determined. Photometric red- 
shifts from SDSS are z — 0.295 ± 0.042 (KD-tree method) and 
z = 0.219 + 0.088 (RF methoc0), both being significantly higher 
than the value (z — 0.045) we obtained directly from spectra 
of LSQ13cuw. Given the unreliability of photometric redshifts, 
combined with the fact that the spectrum is reminiscent of an old 
population, and therefore unlikely to give rise to core-collapse 
SNe, we deem it improbable that SDSS J023958.22-083123.5 
hosted LSQ13cuw. 

In the SDSS-DR10 database, there is a faint extended source, 
SDSS J023957.37-083123.8, at the SN location. Its photomet¬ 
ric redshifts are 0.014 ± 0.082 (KD-tree method) or 0.049 ± 
0.028 (RF method). The latter would match the redshift of the 
SN within the errors. However, the SDSS DR10 r'-band magni¬ 
tude of SDSS J023957.37-083123.8 (m r . = 20.88 ± 0.27) is not 

6 A description of the KD-tree and the RF method can be found at the 
following url: https://www.sdss3.org/dr8/algorithms/photo-z.php. 
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Fig. 1: LSQ13cuw and its environment. Short dashes mark the 
location of the supernova at a J 2000 = 02 /! 39”'57'.35, tfnooo = 
-08°31'24".2. The nearest galaxy, marked “a”, originally pre¬ 
sumed to be the host, is SDSS J023958.22-083123.5. The num¬ 
bers mark the positions of the sequence stars (see also Table lB. it 
used for the photometric calibrations, r'-band image taken on 
2013 November 29, 31.6 days after explosion. 



Aobs(A) 


Fig. 3: Likely host galaxy of LSQ13cuw. A faint extended source 
was detected within 0".5 of the reported position of SDSS 
J023957.37-083123.8,1".3 NEfromLSQ13cuw. The source ap¬ 
pears elongated in the N-S direction with two bright regions at 
each end that are marked with “b” and “c”. “a” marks the po¬ 
sition of the galaxy SDSS J023958.22-083123.5 that was orig¬ 
inally presumed to be the host. Deep r'-band imaging obtained 
on 2014 September 24, 316.6 days after explosion (rest frame). 


ages, or are components of a single irregular galaxy). The mag¬ 
nitude was measured using aperture photometry, and an aperture 
of 3" was used which encompasses the entire source (i.e. regions 
“b” and “c”). The magnitudes of the source are m g > = 22.38 ± 
0.13, m r ' = 22.38 ± 0.13 and m? = 22.37 ± 0.32. If the bright 
region on the north side of the source (“c” in Figure [3} is ex¬ 
cluded, then the r' and /'-band magnitudes are fainter by ~ 0.6 
mags. The source is unresolved in the gr'-band. Given the spatial 
coincidence, we favour SDSS J023957.37-083123.8 as the host 
galaxy of LSQ13cuw, with M r - = -14.06 ± 0.18 (for the entire 
source “b” and “c”), derived using our redshift estimate above. 


Fig. 2: Spectrum of the galaxy SDSS J023958.22-083123.5 des¬ 
ignated with “a” in Figures[T|and[3] 


consistent with the limiting magnitude of the LSQ pre-explosion 
images (nv > 21.5). Indeed it is flagged as being unreliable, so 
any agreement between the SDSS photometric redshift and our 
own spectroscopically-determined redshift is likely to be purely 
coincidental. In order to reliably determine the magnitude of this 
source, we obtained deep imaging with the NTT (gr#782, r#784 
and ;#705 filters) in September and October, 2014. We took 6 
x 200 s exposures in each band, and reduced them using the 
PESSTO pipeline, as detailed in Section I2T1 After aligning and 
combining the images, we found that the SN had faded beyond 
detection, and a faint extended source was detected within 0".5 
of the reported position of SDSS J023957.37-083123.8,1" .3 NE 
from LSQ13cuw. The source appears elongated in the N-S direc¬ 
tion with two bright regions at each end (marked with “b” and 
“c” in Figure [3j it is unclear whether the two bright regions are 
two separate galaxies which are marginally resolved in the im- 


2.3. Photometry 

Table Q] shows the log of imaging observations. The light curve 
is presented in Figure[4] 

The fact that LSQ13cuw was discovered only two days after 
the last non-detection allows us to put tight constraints on the 
explosion epoch. Using these and a low-order polynomial fit to 
the pre-peak photometry we derive the date of explosion to be 
MJD 56593.42 ± 0.68. This makes LSQ13cuw one of the earliest 
discovered SNe II-L. 

LSQ 13cuw rises to peak in 15.9+ 1.2 d (rest frame) with the 
fitted epoch of maximum light in the r'-band being MJD 56610.0 
±1.0 and reaching a magnitude of 18.4 ± 0.1. The r'-band light 
curve then declines linearly until about 26.4 d (rest frame) after 
the peak magnitude, then the decline rate slightly slows down. A 
similar decline to that in the r'-band was also observed in the g'- 
and /'-bands. 

As previously mentioned, the photometric and spectroscopic 
classifications of SNe are not mutually exclusive. However, the 
focus of this work is on explosion-energy driven SNe, so we 
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Table 1: Photometry of LSQ13cuw 


Date 

MJD 

Epoch* 

rest-frame 

mag 

V 

error 

K-corr** 

mag 

9’ 

error 

K-corr** 

mag 

r’ 

error 

K-corr** 

mag 

i' 

error 

K-corr** 

Telescope 

+Instrument 

2013 Oct 26 

56591.31 

-2.02 

- 

- 

- 

- 

- 

- 

>21.19 

- 

- 

- 

- 

- 

LSQ ST 

2013 Oct 28 

56593.31 

-0.11 

- 

- 

- 

- 

- 

- 

>21.21 

- 

- 

- 

- 

- 

LSQ ST 

2013 Oct 30 

56595.13 

1.64 

- 

- 

- 

- 

- 

- 

20.57 

0.36 

<0.03 

- 

- 

- 

LSQ ST 

2013 Oct 30 

56595.21 

1.71 

- 

- 

- 

- 

- 

- 

20.34 

0.24 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 01 

56597.44 

3.85 

- 

- 

- 

- 

- 

- 

19.53 

0.13 

- 

- 

- 

- 

CRTS CSSST 

2013 Nov 03 

56599.19 

5.52 

- 

- 

- 

- 

- 

- 

19.59 

0.23 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 05 

56601.11 

7.36 

- 

- 

- 

- 

- 

- 

19.16 

0.14 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 05 

56601.19 

7.43 

- 

- 

- 

- 

- 

- 

19.15 

0.13 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 07 

56603.11 

9.27 

- 

- 

- 

- 

- 

- 

18.76 

0.11 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 07 

56603.19 

9.35 

- 

- 

- 

- 

- 

- 

18.84 

0.13 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 09 

56605.10 

11.17 

- 

- 

- 

- 

- 

- 

18.55 

0.13 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 09 

56605.18 

11.25 

- 

- 

- 

- 

- 

- 

18.63 

0.11 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 10 

56606.79 

12.79 

- 

- 

- 

- 

- 

- 

18.64 

0.03 

- 

- 

- 

- 

CRTS CSSST 

2013 Nov 10 

56606.80 

12.80 

- 

- 

- 

- 

- 

- 

18.54 

0.03 

- 

- 

- 

- 

CRTS CSSST 

2013 Nov 10 

56606.80 

12.80 

- 

- 

- 

- 

- 

- 

18.80 

0.03 

- 

- 

- 

- 

CRTS CSSST 

2013 Nov 10 

56606.81 

12.81 

- 

- 

- 

- 

- 

- 

18.46 

0.03 

- 

- 

- 

- 

CRTS CSSST 

2013 Nov 13 

56609.09 

14.99 

- 

- 

- 

- 

- 

- 

18.41 

0.17 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 13 

56609.17 

15.07 

- 

- 

- 

- 

- 

- 

18.41 

0.17 

<0.03 

- 

- 

- 

LSQ ST 

2013 Nov 24 

56620.07 

25.50 

19.44 

0.05 

-0.02 

- 

- 

- 

- 

- 

- 

- 

- 

- 

NTT+EFOSC2 

2013 Nov 25 

56621.21 

26.59 

- 

- 

- 

19.35 

0.09 

0.01 

19.17 

0.06 

0.03 

19.05 

0.08 

-0.08 

LCOGT 

2013 Nov 26 

56622.19 

27.52 

- 

- 

- 

- 

- 

- 

19.34 

0.09 

- 

- 

- 

- 

NTT+EFOSC2 

2013 Nov 26 

56622.28 

27.61 

- 

- 

- 

19.46 

0.12 

0.01 

19.33 

0.10 

0.04 

19.13 

0.13 

-0.08 

LCOGT 

2013 Nov 28 

56624.12 

29.37 

- 

- 

- 

19.63 

0.07 

0.02 

19.32 

0.05 

0.06 

19.22 

0.06 

-0.09 

LCOGT 

2013 Nov 29 

56625.03 

30.24 

- 

- 

- 

19.67 

0.09 

0.02 

19.40 

0.03 

0.06 

19.23 

0.04 

-0.10 

LT+IO:0 

2013 Dec 01 

56627.19 

32.31 

- 

- 

- 

- 

- 

- 

19.62 

0.08 

- 

- 

- 

- 

NTT+EFOSC2 

2013 Dec 07 

56633.90 

38.73 

- 

- 

- 

20.40 

0.10 

0.03 

19.95 

0.15 

0.13 

19.74 

0.09 

-0.15 

LT+IO:0 

2013 Dec 16 

56642.86 

47.30 

- 

- 

- 

20.84 

0.36 

0.04 

20.27 

0.18 

0.20 

20.11 

0.13 

-0.22 

LT+IO:0 

2013 Dec 21 

56647.85 

52.07 

- 

- 

- 

21.02 

0.12 

0.04 

20.37 

0.09 

0.24 

20.43 

0.11 

-0.26 

LT+IO:0 

2013 Dec 22 

56648.12 

52.33 

20.99 

0.07 

-0.04 

- 

- 

- 

- 

- 

- 

- 

- 

- 

NTT+EFOSC2 

2013 Dec 23 

56649.94 

54.07 

- 

- 

- 

20.84 

0.25 

0.04 

20.47 

0.49 

0.26 

20.20 

0.26 

-0.28 

LT+IO:0 

2013 Dec 24 

56650.84 

54.93 

- 

- 

- 

20.89 

0.13 

0.04 

20.54 

0.12 

0.27 

20.38 

0.12 

-0.29 

LT+IO:0 

2013 Dec 25 

56651.83 

55.88 

- 

- 

- 

20.91 

0.14 

0.04 

20.70 

0.15 

0.27 

20.49 

0.13 

-0.30 

LT+IO:0 

2013 Dec 26 

56652.84 

56.84 

- 

- 

- 

21.07 

0.11 

0.05 

20.59 

0.07 

0.28 

20.43 

0.09 

-0.31 

LT+IO:0 

2014 Jan 11 

56668.98 

72.29 

- 

- 

- 

21.53 

0.49 

0.06 

20.97 

0.39 

0.41 

20.66 

0.29 

-0.50 

LT+IO:0 

2014 Jan 13 

56670.98 

74.20 

- 

- 

- 

- 

- 

- 

20.81 

0.28 

0.43 

20.71 

0.28 

-0.53 

LT+IO:0 

2014 Jan 24 

56681.08 

83.86 

22.04 

0.22 

-0.07 

- 

- 

- 

- 

- 

- 

- 

- 

- 

NTT+EFOSC2 

2014 Jan 30 

56687.08 

89.60 

22.22 

0.35 

-0.07 

22.06 

0.13 

- 

21.13 

0.10 

- 

21.45 

0.17 

- 

NTT+EFOSC2 

2014 Feb 08 

56696.85 

98.95 

- 

- 

- 

- 

- 

- 

21.92 

0.15 

0.63 

21.77 

0.18 

-0.95 

LT+IO:0 


*The phase is given relative to our estimate of the epoch of explosion on MJD 56593.42. Limiting magnitudes are to 3 cr. **K-corrections were calculated by first using the 
LSQ13cuw spectroscopy and the snake code (SuperNova Algorithm for K-correction Evaluation) within the S3 package (Inserra et al. in prep.). These values were then 
interpolated to all epochs of photometric observations. The K-corrections are not included in the photometry, as up to an epoch of about 30 days after explosion the derived 
K-corrections are smaller or similar to the error in the photometry whereas our analysis focusses primarily on the early epochs. LSQ = La Silla-QUEST; ST = ESO 1.0m 
Schmidt Telescope; CRTS = Catalina Real-time Transient Survey; CSSST = Catalina Sky Survey Schmidt telescope; LT = Liverpool Telescope; NTT = New Technology 
Telescope; LCOGT = Las Cumbres Observatory Global Telescope Network lm telescope. We note that the NTT filters are not standard Sloan filters and that CRTS data are 
unfiltered, see also Section [2~T1 The transmission curves for the respective filters are available at: LT g'r'i': http://telescope.livjm.ac.Uk/TelInst/Inst/IOO/i LCOGT g'r'i 
http://lcogt.n et/configdb/filter/; NTT E#641, g# 782, r#784, /#705: http://www.eso.org/sci/facilities/lasilla/instruments/efosc/inst/Efosc2Filters.html; LSQ wideband filter: 
iBaltav et aD (120131) . Figure 2. 
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Days since explosion (rest frame) 

0 20 40 60 80 100 



MJD 


Fig. 4: Vg'r'i' light curves of LSQ13cuw. The arrows represent 
pre-discovery limits. The vertical lines indicate epochs of spec¬ 
troscopy. We note that the NTT filters are not standard Sloan 
filters and that CRTS data are unfiltered (Section [2~T1) . 



Arest(A) 


wish to exclude any objects of type Ibc/IIb/IIn even if they have 
a linearly declining light curve. Of these, the type lib SNe are 
probably the most problematic. Some SNe lib show a_distinct 
primary peak, but most do n ot (e.g. SN 20 08ax, Pastore lo et all 
12008k iTsvetkov et alll2009l or SN 201 ldh, lErgon et al.ll2014bl) . 


As LSQ13cuw was discovered early, it is unlikely that a primary 
peak was missed. We defer a discussion of relative durations of 
rise times to Section l3~4l 

During the post-peak decline phase, the /(-band light curves 
of SNe lib show a similar behaviour to SNe II-L in that they have 
a relatively linear decline. However, a characteristic property of 
SNe lib is a very distinct flattening in the blue bands ( U,B,V,g ; 
see e.g. Figure 2 in lErgon et al.ll2014ah . 

The flattening around 40 days post-explosion in Type libs is 
related to the transition from the 56 Ni-powered diffusion phase 
to the 56 Ni-powered tail (quasi-nebular) phase. In a II-L or II-P 
SN the light curve in the first ~ 100 days is instead driven by the 
diffusion of explosion-deposited energy. 


In order to determine whether LSQ13cuw bore more of 
a resemblance to a type lib or II-L SN, we measu red the 


decline rates of the Type lib SNe 1993J (iRiperoetal. 

1993; 

Corwin et al.l 19931; Centurion et al.l 1993| Hanzl et al. 

1993; 

Tweedy et al. 1993; Dumont et al. 1993? Pressberger et al 

1993; 

Lewis et al.; i^fMclBarbon et al.lll995). 2008ax, (Pastorello et al. 

2008; Tsvetkov et alj 20091 and 201 ldh (jErgonetal.il 2 

oh™. 

of the Type II-L SNe 1979C ([Balinskaia et al. 

19S(jt 

de Vaucouleurs et al.l 1981; Barbon et al. 1982h, 

1990K 

(Cappellaro et al.l 1995b and 2001fa (Faran et al. 2014bl) as well 


as of LSQ13cuw. For the lib SNe, we measure initial post-peak 
V-band decline rates of 6.4-7.3mag/100d, which then flatten 
significantly to declines rates of 1.4-2.2 mag/100 dQ 


SN 1979C declines with at a rate of 3.6 ± 0.2 mag/100d. 
As it was only discovered relatively late it is not clear whether 
it also had a somewhat steeper initial decline from maximum. 


7 Within the errors, our measurements are consistent w ith the tail- 
decline rates of 1.8-1.9 mag/100 d in the V-band reported bv lErgon et all 
li2014ab for the three SNe Type lib 1993J, 2008ax, and 201 ldh. 


Fig. 5: LSQ13cuw spectroscopy. *Flux normalized to the maxi¬ 
mum Htr flux and smoothed for better visibility of the features. 
The exact normalizations and smoothing factors are: flux/6.9 
smoothed by a factor of 3 for the +25 d spectrum; flux/4.3 
smoothed by a factor of 5 for +25 d; flux/4.3 not smoothed for 
+32 d; flux/4.0 smoothed by a factor of 3 for +52 d; flux/2.5 
smoothed by a factor of 3 for +84 d. 


SNe 1990K and 2001 fa both display a slight shoulder of about 
20-30 days in their decline, but then fall even more steeply. SN 
1990K has an average decline rate of 4.1 +0.1 mag/100 d. SN 
2001fa initially declines at a rate of 5.0 + 0.4mag/100d but then 
flattens to an over all decline rate of 4.1 + 0.1 mag/100 d. 

In contrast, LSQ13cuw initially declines by 5.8 + 
0.6 mag/100d the g '-band and flattens to a decline rate of 3.0 + 
0.5 mag/100 d in the tail. Over-all the g'-band decline behaviour 
of LSQ13cuw seems to match the decline rates measured for 
the Type II-L SNe and we conclude that our classification of 
LSQ13cuw as a Type II-L SN is justified. 

We are aware that both SNe lib and II-L show a broader 
range in light curve decline properties than is outlined here, and 
other studies previously mentioned above have focussed on the 
analysis of Type II decline rates. 

2.4. Spectroscopy 

Table [2] shows the journal of spectroscopic observations. The 
fully reduced and calibrated spectra of LSQ13cuw are presented 
in Figure [5] They are corrected for reddening (E(B - V) = 
0.023 mag) and redshift (z = 0.0453). 

The strongest feature in the spectrum is Her, which would be 
matched by the corresponding feature at ~ 4850 A in the blue to 
be H/3. 

In the earliest spectra, both the Her and H/j profiles are rel¬ 
atively weak but they become stronger at later epochs (see also 
Table [3]) and sho w at best o nly a very weak absorption com¬ 
ponent. Recently iGutierrez et al.l (12014b analysed a sample of 
hydrogen-rich Type II SNe and found that the ratio of absorp- 
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Table 2: Journal of spectroscopic observations 


Date 

MJD 

Epoch* 

rest-frame 

Wavelength 
range in A 

Telescope 

+Instrument 

Resolution 

A 

2013 Nov 24 

56620.08 

+25.49 

3550-8670 

NTT+EFOSC2+Grl3 

18.0 

2013 Nov 26 

56622.21 

+27.53 

3600-7150 

NTT+EFOSC2+Grl 1 

13.8 

2013 Dec 01 

56627.22 

+32.32 

3600-8060 

NTT+EFOSC2+Grl 1,16 

14.0 

2013 Dec 22 

56648.15 

+52.33 

3600-8860 

NTT+EFOSC2+Grl 1,16 

14.0 

2014 Jan 24 

56681.13 

+83.86 

3800-8380 

NTT+EFOSC2+Grl3 

17.7 


*Rest frame epochs (assuming a redshift of 0.0453) with respect to the assumed explosion date of 56593.42 (MJD). The resolution 
was determined from the FWHM of the sky lines. 


tion to emission in the Ha P-Cygni profile is smaller for brighter 
and faster declining light curves and for SNe with higher Ha 
velocities. They attribute this to the mass and density profile of 
the hydrogen envelope. In this picture we would expect SNe II- 
L to have smaller Ha absorptions than SNe II-P. The fact that 
LSQ13cuw agrees with these results is en cour aging. The sample 
of SNe II-L presented by iFaran et aD ( 2014bl) is also in accord 
with this conclusion. 

The shape of the Ha line profile evolves quite drastically in 
most Type II SNe. In the earliest spectrum of LSQ13cuw Ha is 
weak, however a few days later the Ha line becomes stronger 
and has a fla t top. This might be a p ossible indication of a de¬ 
tached layer (IJefferv & Branch!! 1990 ). By +52 d, the seemingly 
flat top of the Ha feature seen at +32 d has vanished, leaving a 
sawtooth shaped profile in its wake. The peak of the Ha line is 
blue-shifted by ~3600 km s -1 compared to the rest wavelength. 
Similarly, the Ha line in the +84 d spectrum is blue shifted by 
~2000kms _1 . 

A blue-shifted peak in the Ha profile was already noted 
for SN 1979C. IChugail d 19851) suggest this might be due to an 
opaque core screening the atmospheric layers in the expand¬ 
ing envelope in which the Ha line is formed. [Anderson et al.l 
( 2014a) have investigated the occurrence of blue shifted emis¬ 
sion peaks in 95 Type II supernovae and come to the conclu¬ 
sion that it is a generic property, observed for the majority of 
objects in their sample. They perform non-LTE time-dependent 
radiative-transfer simulations and show that the shape of the Ha 
line profile is a result of the steep density profile of the H layers 
of the ejecta. Since the opacity is dominated by electron scat¬ 
tering and also the line emissions mainly come from the region 
below the continuum photosphere the probability of observing 
blue-shifted line photons is higher. This effect decreases with 
time as the ejecta expand and the density decreases. 

We can also discern a feature at about 5900 A. This is typi¬ 
cally attributed to a blend of He 1/15876 and Na i D 4/15890,5896. 
As with the hydrogen line profiles, it only shows a very weak 
absorption. At about 7100 A there seems to be another emission 
feature, which shows weakly in the +32 d spectrum and a bit 
stronger in the +52 d spectrum. We suspect this to result from 
He i 47065. 

In Figure[6]we compare the spectra of LSQ13cuw with spec¬ 
tra from other Type II SNe. For all SNe the epoch is given rel¬ 
ative to the estimated date of explosi on. The exception is SN 
1980K (e.g. lUomoto & Kirshnerll 19861) for which no constraints 
on the explosion epoch are available and the epochs are therefore 
given relative to the first detection. 

In the earliest spectrum of LSQ13cuw we see a few weak 
features similar to the early spectrum of the Type II-L SN 1980K. 
The flat top of the Ha line of LSQ13cuw a few days later has 
previously also been observed in the early spectra of the Type 



^resttA) 


Fig. 6: LSQ13cuw spectra in compariso n with 
the SNe 1980K (II-L; 

(lib; J c tfcrv ct'al.i 19941 


Uomoto & Kirshnerl 


Barbon et a ] 1995 


198 


ectra of 
1993J 


Fransson et al.l 


20051), 2009hd (II- L; lElias-Rosa et al.l 1201 ll) and 201 ldh (lib; 


Ergon et aTl l2014bh . The epochs are given relative to the esti 


mated date of explosion, with the exception of SN 1980K for 
which the epochs are given relative to the first detection. The 
dotted lines correspond to the positions of Hei 45015, 45876, 
46678 and 47065. *Flux normalized to the maximum Ha flux 
and smoothed for better visibility of the features. 


lib SN 1993J. However, in the further evolution of SN 1 993J the 
Ha emission splits into two components. iMatheson et akl (2 00(1 ) 
attributed this to clumping in the ejecta of SN 1993J. We see no 
evidence of this in the spectra of LSQ13cuw0 

The shape of the Ha profile is quite similar in the +25 and 
+32 d spectra of LSQ13cuw, the +19 d spectrum of SN 1993J, 

8 The small absorption in the centre of the Ho- profile of LSQ13cuw is 
a telluric feature. 
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Table 3: Emission equivalent width measurements for LSQ13cuw. 


Date 

MJD 

Epoch* 

rest-frame 

Her 

EWin A 

H/J 

EW in A 

Hei/NaiD“ 
EW in A 

2013 Nov 24 

56620.08 

+25.49 

-41 + 1 

-16 + 7 

-34 + 4 

2013 Nov 26 

56622.21 

+27.53 

-109+10 

-18 + 2 

-41+5 

2013 Dec 01 

56627.22 

+32.32 

-152 + 24 

-57 + 8 

-47 + 8 

2013 Dec 22 

56648.15 

+52.33 

-480 + 73 

-170 + 37 

-146 + 27 

2014 Jan 24 

56681.13 

+83.86 

-800+ 190 

-242 + 85 

-276 + 50 


‘Rest frame epochs (assuming a redshift of 0.0453) with respect to the assumed explosion date of 56593.42 (MJD). “He i 45876 
and Na i D 445890,5896 blend. 


and the +19 d spectrum of 2009hd, even though they represent 
epochs that are two weeks apart. This might indicate that differ¬ 
ent SNe go through the different evolutionary stages on different 
time scales. 

The sawtooth shape in the +52 d spectrum the Her feature 
in LSQ13cuw strongly resembles the Her profile of the >55 d 
spectrum of SN 1980K, which is similarly asymmetric. 

Compared to the other Type II-L SNe like 1980K and 
2009hd, it seems that in LSQ13cuw the Her feature is generally 
broader. Only the SN 1993J spectra display similarly broad Her 
profiles. This might be an indication of a somewhat broader ve¬ 
locity distribution. 

Apart from Her we also see H/3 in the various SNe as well as 
the blend of He 1 45876 and Na i D 445890,5896. The He 1 47065 
line does not show in any of the comparison spectra. The feature 
between 7000 and 7300 A in the late spectrum of SN 1993J is 
probably due to the [Ca n] 447291,7324 doublet rather than He i. 

Given the presence of a feature near 45880 A, close to Hei 
45876 in the 25 and 32 d spectra of LSQ13cuw (Figures 5, 6), 
we must consider whether a spectroscopic classification of Type 
lib is warranted. A closer inspection reveals marked differences 
to Type lib spectra (Figure 6). We note that prototypical lib SNe 
such as SNe 1993J and 201 ldh show optical signatures of helium 
about two weeks post-explosion (e.g. ISahu et al.ll2013l) with hy¬ 
drogen lines becoming progressively weaker with time. In con¬ 
trast, LSQ13cuw maintains a strong feature due to Her at epochs 
as late as 84 d, while the putative 45876 feature is no longer 
discernible. Although we cannot rule out the emergence of He i 
46678 in the red wing of Her giving rise to the pronounced asym¬ 
metry, the relative strengths of the hydrogen versus helium lines 
argues against a firm Type lib classification. In Section 3.4 we 
show that the rise time behaviour of bona fide lib SNe is also 
markedly different from Type IIP-L SNe. 

Finally, it is interesting to note that compared to the other 
SNe the spectra of LSQ13cuw show almost no, or only ex¬ 
tremely weak lines of intermediate mass elements between 4000 
and 5500 A, which we would expect to appear a few weeks after 
explosion. 


3. Results and discussion 

3.1. The sample 


The scarcity of SNe II-L poses a challenge in the quest to com¬ 
pare LSQ13cuw to other SNe II-L. There are a few well-known 
objects in this class l ike SN 1979C (e.g.lBalinskaia et al. 1980t 


Panagia e t al.l 1980l: 


r 


1981 Barbon et al 


de Vaucouleurs et ah 

119821) . SN 1980K 


198 


(e.g. 


Thompson! 19821), SN 199 0K (ICannellaro et al 


.B ranch et al. 

Butal Il982blat 


19951) and SN 


1996L ( Benetti et ahll 19991) . however while we can use their data 


to compare the post-maximum behaviour, none of the aforemen¬ 
tioned objects was observed during its rise to peak. 


iFaran et ah (2014b) recently published a set of eleven Type 
II-L SNe, three of which have good pre-maximum photometry 
and a well constrained explosion epoch, that we can therefore 
add to our sample: SNe 2001 fa, 2003hf and 2005dq. 

For SNe II-P the situation is considerably better with quite 
a few SNe II-P having bee n discovered and ob s erved at early 
stages: SN 1999 em (e . g, Hamuy et all 1200 it [Leonard et all 
l2002at iLeonardl 12002: iFlmhamdi et al.l 120031). SN 2 004et 


(e.g. Li et al. 2005T Sahu et al. 20Qg 


SN 2 005cs (e.g. Pastorello et alT 2006fc iTsvetkov et al.l 120061 : 

iBrown et al.ll2007l) . to name a few, but see also Table[4] 


Kotak et al 


and 


Our analysis is mainly performed in the r'/R- band, where 
all objects in our sample have reasonably good quality data. In 
the following, we will refer to the r'/R -band for all observables 
unless explicitly stated otherwise. 

Our sample and the light curve properties are presented in 
Table|4] Figure[7]shows a comparison of pre- and post-peak light 
curves for selected SNe from our sample. 


3.2. Rise time parametrization 


We investigate the rise times of 20 Type II-P and II-L SNe with 
explosion epochs constrained between + 0.3 days and + 8 days. 


For most SNe the epoch of explosion was adopted from the 
literature. In the cases where no estimate was given we derived 
the explosion epoch by taking the average MJD between the first 
detection and the last pre-discovery non-detection (SNe 1979C, 
2000dc, 200Icy, 1999gn, 1999gi, 2009hd, 2001do, 2001fa). For 
SNe 2004du and 2004et no pre-discovery non-detections were 
reported, we therefore used a low-order polynomial fit to the 
observed data to estimate the explosion epoch. Finally for SN 
1999em multiple estimates for the explosion epoch using the 
expanding photosphere method are available. For our analy- 
si s we used the mean value dervied from t h e best estimate s 
of lHamuv et alJ (l2 001l).lLeonard et all (l2002al) :lL eonar dj (i200l . 
and Elmhamdi et alJ ( 20031) . 


In order to derive the rise time in an internally self-consistent 
manner, we define an epoch termed “end-of-rise” which is 
the epoch at which the /-'-band magnitude rises by less than 
0.01 mag. d -1 . This was estimated by fitting a low-order polyno¬ 
mial to the data, with an iteratively chosen step-size in time. We 
found that this approach allowed the inclusion of SNe with dif¬ 
ferent light curve morphologies (Fig. [Tat : while the light curves 
of some SNe (such as LSQ13cuw) have a clear peak and a well 
defined maximum, others (such as SN 1999em) display no clear 
peak before settling onto the plateau or (e.g. SN 2005cs) have a 
rising /-band light curve after an initial weak maximum, which 
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Table 4: Comparison of Type II SN properties 


SN 

Type* 

Host galaxy 

E(B 

Galactic 

- V) 

Host 

/J** 

t Q (MJD)“* 

End of rise 
absolute magnitude 

Rise time in 
days (rest frame) 


SN 1990E 

II-P 

NGC 1035 

0.022 

0.46 + 0.10 

31.48 + 0.22 

47934.56 ± 2.56 

-17.69 + 0.34 

9.6 ± 2.6 

a,b,c,d,e 

SN 1999em 

II-P 

NGC 1637 

0.04 

0.06 

30.23 + 0.14 

51475.87 + 2.59 

-16.90 + 0.14 

8.3 ±2.7 

a,f,g,h,i 

SN 1999gi 

II-P 

NGC 3184 

0.017 

0.19 + 0.09 

30.54 ± 0.29 

51514.84 + 0.34 

-16.75 + 0.38 

11.5 + 0.6 

a,b,j,k,l 

SN 2002gd 

II-P 

NGC 7537 

0.059 

- 

32.87 + 0.35 

52551.5 ±2.0 

-16.10 + 0.35 

8.1 ±2.2 

a,m 

SN 2004du 

II-P 

UGC 11683 

0.084 

- 

33.94 + 0.13 

53225.6 ± 1.1 

-17.67 + 0.13 

10.9+ 1.5 

a,b,n 

SN 2004et 

II-P 

NGC 6946 

0.34 

0.07 ± 0.07 

28.37 + 0.15 

53270.61 ± 0.60 

-17.15 + 0.24 

8.6 ± 0.9 

a,b,o,p,q,r 

SN 2005cs 

II-P 

NGC 5194 

0.032 

0.08 ± 0.04 

29.26 ± 0.33 

53548.5 ±0.5 

-15.12 + 0.35 

3.1 ±0.5 

a,s,t,u,v,w,x 

SN 2006bp 

II-P 

NGC 3953 

0.027 

- 

31.33 + 0.19 

53833.7 ±0.3 

-16.46 + 0.19 

6.7 ± 0.4 

a,b,w,y 

SN 2010id 

II-P 

NGC 7483 

0.054 

- 

33.12 + 0.48 

55452.3 ±2.0 

-14.37 + 0.48 

6.2 ± 2.0 

a,b,z 

SN 2012A 

II-P 

NGC 3239 

0.029 

- 

29.96 + 0.15 

55933.0 

-16.33 + 0.15 

8.1 ± 3.0 

a,A 

SN 2000dc 

II-P/II-L 

ESO 527-G019 

0.071 

- 

32.93 + 0.14 

51761.8 ±3.5 

-17.25 + 0.19 

10.5 + 4.9 

a,b,n,B 

SN 200Icy 

II-P/II-L 

UGC 11927 

0.185 

- 

33.01 + 0.12 

52085.0 ±5.5 

-17.65 + 0.12 

5.4 ±5.5 

a,b,n,B,C 

SN 200 ldo 

II-P/II-L 

UGC 11459 

0.170 

- 

32.35 + 0.15 

52133.2 ±2.0 

-17.27 + 0.15 

8.1 ±2.0 

a,b,n,B 

SN 2013ej 

II-P/II-L 

NGC 0628 

0.061 

- 

29.77 ± 0.28 

56496.95 ± 0.92 

-17.60 + 0.28 

12.9+ 1.4 

a,b,D,E,F,G,H,I 

SN 1979C 

II-L 

NGC 4321 

0.023 

0.16 + 0.04 

30.30 + 0.15 

43970 ± 8 

-18.72 + 0.37 

23.9 + 8.0 

a,b,J,K,L 

SN 2001 fa 

II-L 

NGC 0673 

0.069 

- 

34.09 ± 0.27 

52197.9 ±2.5 

-18.13 + 0.27 

7.8 ±2.6 

a,b,B,M 

SN 2003hf 

II-L 

UGC 10586 

0.019 

- 

35.55 ± 0.04 

52863 ± 2 

-18.46 + 0.05 

13.2 + 2.1 

a,b,B 

SN 2005dq 

II-L 

UGC 12177 

0.070 

- 

34.77 ± 0.03 

53608 ± 4 

-17.50 + 0.03 

9.1 ± 3.9 

a,b,B 

SN 2009hd 

II-L 

NGC 3627 

0.030 

1.20 

29.86 ± 0.08 

55001.6 ±3.9 

-17.11 + 0.16 

15.9 + 4.1 

a,b,N,O.P 

LSQ13cuw 

II-L 

SDSS J023957.37-083123.8 

0.023 

- 

36.38 + 0.13 

56593.42 ± 0.68 

-18.03 + 0.13 

15.0 ±0.8 

a,b 


‘Type for each SN adopted from literature; **Hq — 72 km s -l Mpc 1 “‘For most SNe the epoch of explosion was adopted from the literature. In the cases where no estimate was 
given we derived the explosion epoch by taking the average MJD between the first detection and the last pre-discovery non-detection (SNe 1979C, 2000dc, 2001cy, 1999gn, 
1999gi, 2009hd, 2001do, 2001fa). For SNe SN 2004du and 2004et no pre-discovery non-detections were reported, we therefore used a low-order poly nomia l fi t to the observed 
dat a to estimate the explosi on epoch. For SN 1999em we adopted the average of the explosion epochs derived by lHamuv et al.1 ( 2001 1), Leonardetah ( 2002al);lLeonar d J2002I) 


a ndlElmh amdi et aU (12003 ) using th e expanding photos ph ere method, a) this paper; b) NASA/I P AC Extragalactic Datab as e: cllPennvna cker et a l.lll990t d) Schmidt et al.ll 1 993b 
el Benetti et al.ll 19941: f) Li 1999c gflHamuv et alJl200l t hi [Leonard et all 2002at Leonardl2 002t if iHm hamdi et al.ll2003fc il Nak ano & Ku shidall 1 9 99t ki Trondaletalj tjj99t 11 


2006; t) Takats & Vinkoll2006l: u) Tsvetkov et alll2006; v) Brown et ap|2007 : ; w) Dessart et all|20f 

>8U) 

Pastorello et al. 

2009; y) Ouimby et al. 2007!: z) Gal-Yam et akEoi 1 

I A) 

Tomasella et al. 2013!; B) Faran et al][2014b; C) Ganeshalinsam et al. 2001 D) Shappee et al. 201 

iEi 

Kim etalj2013 

; F) Dhungana et al. 20131 G) Valenti et al. 2014; H) 



Richmondll2014t I I Bose et al.||2015 ; J) Balinskaia et al.ll 19801: Kl de Vaucouleurs et al.ll 198 it LI Barbon et al.ll 1 982t Ml lFilippenko & Chornockll2001 ; N) Elias-Rosa et al.luOl it 


Ol iKasliwal et al.ll2009b PI lMonardll2009[ 
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Days since explosion (rest frame) 



(a) r'/R -band pre-peak light curves. 


(b) r'/R -band post-peak light curves. 


Fig. 7: r'/R -band light curve comparison of a selection of SNe II-P and II-L from our sample. Magnitudes are plotted relative to the 
“end-of-rise” magnitude (see Section [3~2l i to visualize the range in rise times and decline rates amongst the various SNe. The arrows 
show the limits on the pre-discovery non-detections. The dotted line marks the explosion epoch. 


would lead to an epoch of maximum light well into the plateau 
phase. 

To estimate the error in this process, we linearly extrapo¬ 
lated the fit to the point where the gradient equalled zero. The 
midpoint between this value and either the last measured data 
point, or the “end-of-rise” epoch (whichever one was more dis¬ 
tant) was taken to be a conservative estimate of the error. This 
naturally accounts for how well-sampled the rise to maximum 
or plateau actually was. With the exception of SN 1979C the 
SNe in our sample have a sampling quality of 0.3-1.3 photo¬ 
metric points per day (averaged over ± 3 days around the fitted 
result for the “end-of-rise” epoch), which is reasonably good for 
our purposes. The final error on the rise time is given by adding 
the errors on the explosion epoch and the “end-of-rise” epoch in 
quadrature. Analogously the error on the “end-of-rise” absolute 
magnitude is determined by the errors on the distance modulus, 
the “end-of-rise” magnitude and the total reddening. The error in 
the explosion epoch estimate is combined in quadrature with the 
error in the gradient; the former dominates the combined error 
in all cases. 

Figure [8] shows a comparison of the rise times with the ab¬ 
solute magnitudes at the “end-of-rise” epoch. The filled squares 
represent SNe that were previously classified as II-L, while the 
filled circles were classified as SNe II-P. The empty square and 
circle represent ri se time meas urements for the II-P and II-L 
templates from lFaran et al.l (l2014bl) . respectively. For the abso¬ 
lute magnitu des of these templates we assumed the mean values 
published bv [Richardson et ah ( 2014 3 for the two types, respec¬ 
tively. We note, however, that these are /I-band values. The ar¬ 
row represents a lower limit for the rise time of the Type II-L SN 
2009kr. Unfortunately more than one month passed between the 
last non-detection and the discovery of SN 2009kr dNakano et al.i 
120091 3 and the explosion epoch is therefore not very well con¬ 
strained. Ho wever, pre-peak p h otom etry for SN 2009kr has been 
published bv lElias-Rosa et all (1201 Oh and the difference between 
the discovery and the estimated epoch of “end-of-rise” was used 
as a lower limit for the rise time. We also included the peculiar 
Type II-P SN 1987A and the Type lib SNe 1993J, 2008ax, and 
201 ldh for comparison (see also Table HOI in the Appendix). 



Rise time in days (rest frame) 


Fig. 8: Comparison of “end-of-rise” r'/R -band absolute magni¬ 
tudes and rise times. The filled squares represent SNe that were 
classified as Type II-L, while the filled circles are SNe that were 
classified as Type II-P. The empty square and circle represent 
rise t ime measu rements for the R-band II-L and II-P templates 
from lFaran et al.l (12014hb . respectively. As these templates are in 
a r elative scale, we assume d the mean B-band values published 
by iRichardson et al.l ( 2014 1 as the absolute magnitudes for the 
two SN types, respectively. The arrow represents a lower limit 
for the rise time of SN 2009kr. For the individually labelled SNe 
2005cs, 2009kr and 2012A direct progenitor detections have 
been reported (see also Section l3~5l) . The peculiar Type II-P SN 
1987A and the Type lib SNe 1993J, 2008ax and 201 ldh are in¬ 
cluded for comparison (see Table IB. 2l . 


For these SNe the explosion epochs are constrained to less than 
± 1 d. For the estimate of the rise time the first peak was used in 
the case of SN 1987A and the second peak was used in the case 
of SN 1993J. 
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3.3. Rise time correlations 

An inspection of Table [4] and Fig. [7a] reveals some intrigu¬ 
ing points which we highlight below. First, only considering 
those few objects in our sample that have the most tightly con¬ 
strained explosion epochs (<1 d) and peak magnitudes, we find 

that SN 2005cs, an undisputed type II-P SN, rose to Mr -15 in 

~3 d, while LSQ13cuw, with type II-L characteristics, had a rise 

time of ~15 d and a peak magnitude of M'-18. Objects with 

peak magnitudes fainter than LSQ13cuw (such as SNe 2004et, 
2006bp) have correspondingly shorter rise times. This is borne 
out in Fig. [8] albeit with some scatter. 

Other studies have noted the possibility of such_a trend . Fo r 
examp le, the early discovery of SN 2010id led iGal-Yam et al.l 
( 2011 ) to compare its rise time to that of SNe 2005cs and 
2006bp; they found that SNe 2005cs and 2010id - both faint type 
II-P SNe - had shorter rise times than SN 2006bp. On the other 
hand lFaran et al.l f 2014a!) in a recent study, report no conclusive 
correlation between rise time and peak brightness. 

If the above trend holds over all type II-P, L SNe, then objects 
such as SNe 2000dc and 2001do appear to be true intermediates 
in the sense that they are brighter at peak and have longer rise 
times than the average SN II-P, while they are fainter, and have 
shorter rise times than the average II-L . Intere stingly, these very 
objects were culled bv lPoznanski et al] d2009t) from their sample 
of type II-P SNe to be used as distance indicators owin g to their 
larger-than-average post-peak decline rates, while iFaran et akl 
(120 14bl) classified these as type II-L SNe for their study of type 
II SNe. Another case in point is SN 2013ej : dubbe d a “slow- 
rising” type II-P SN by IValenti et al.l (120141) . while iBose et akl 
020151) favour a “Type II-L” classification based on its post-peak 
decline rate. We find it to be similar to the “fast-declining” SNe 
II-P in that it is relatively bright, and has a relatively long rise 
time compared to other SNe II-P (see Table]]}. 

In Fig. [8] we show the location of the well-known, but pe¬ 
culiar type II SN 1987A as a reference point. It is faint at peak 
relative to the other objects, but has a rise time of ~10d to the 
first peak when estimated in a manner identical to the rest of our 
sample. In the light of our earlier discussion on the possibility 
of lib SNe showing linearly declining light curves, we include 
three well-observed type lib SNe in Fig. [8] With 22.3 ± 0.9, 21.6 
± 0.3 and 20.0 ± 0.5 days the Type lib SNe 1993J, 2008ax and 
201 ldh have a significantly longer rise times than the SNe in 

our sample at a similar brightness (-17.4 to -17.9 mag). The 

only SN with a similar rise time is SN 1979C; however, since 
its explosion epoch is not very well constrained and its rise time 
error therefore correspondingly larger (± 8 days). We note that 
the rise time of SN 1993J, which has a similar luminosity as 
LSQ13cuw is almost one week longer than that of LSQ13cuw 
(15.0 + 0.8 days). 

We next consider whether the type II-P SNe rise times are 
different from those of type II-L SNe. Taken at face value, the 
weighted average of each subsample seems to indicate that there 
is a difference: 7.0 + 0.3 d (II-P) versus 13.3 ± 0.6 d (II-L). How¬ 
ever, we note that the most extreme objects in terms of their rise 
time duration, also have the most tightly constrained explosion 
epochs. We might therefore expect that when objects (of either 
formal subtype) of intermediate peak brightness become avail¬ 
able, a similar treatment is likely to result in an overlap between 
average rise times. We return to this point below. For compari¬ 
son, the rise time measurements for th e II-P and II-L light curves 
templates constructed from the IFaran et ahl (1201 4bl) sample us¬ 
ing the same method, we find rise times of 8.9 + 0.5 d (II-P) and 
11.9 + 0.5 d (II-L). These are shown as open symbols in Fig¬ 


ure [8] Both samples are small and have objects in common, but 
there are significant differences in the make up. For instance, the 
inclusion of SNe known to be int eract ing with a dense circum- 
stellar medium (e.g. SN 2008fq. iTaddia et afl[2013 ). but which 
have linearly declining light curves, could arguably be included 
in studies that consider the post-peak properties of type II-L SNe, 
but are likely to bias rise time estimates in ways not considered 
here. 

There has been much discussion in the literature regarding 
the bimodality or o therwise in the properties of type II-P versu s 
type II-L SNe (e.g. lAnderson et al.ll2014btlSanders et al.ll2014i) . 
Motivated by our findings above, we suggest that an alternative 
way of subdividing these two groups of SNe would be to con¬ 
sider their rise time / peak brightness values if available: SNe 
with longer rise times will generally be brighter at peak. A corol¬ 
lary would be that they would also decline faster post-peak; i.e., 
these would be the undisputed type II-L SNe in the traditional 
classification. Currently, constraints on explosion epochs to bet¬ 
ter than 1 d are not routinely available. However, this is a tech¬ 
nological shortcoming that future planned facilities will remove. 

3.4. Rise time modelling and comparison with observations 

In this subsection we attempt to put our findings into context 
with theoretical expectations. We have restricted ourselves to fo¬ 
cussing on the main results here, with full details of the models 
and associated assumptions presented in the Appendix. 

The optical light curve can be thought of as the bolometric 
light curve multiplied by the fraction of the bolometric light that 
emerges in the optical. As this fraction usually is a faster evolv¬ 
ing function than the bolometric luminosity, the optical peak will 
be close to the peak of the fraction function. We are therefore 
justified in comparing to r'- or R-band observations. 

We performe d a simple pa r amete r stu dy using the analyti- 
cal jnodels from lArnetti ( 198(1 Il982l) and iRabinak & Wax in aril 
(1201 ll) (see Section [AT} assuming a blackbody spectral energy 
distribution and no energy input from 56 Ni in the ejecta. The 
early bolometric light curves were calculated while varying the 
progenitor radius between 100 and 900 R©, the explosion energy 
between 0.1 xlO 51 and 5xl0 51 erg and the ejecta mass between 3 
and 20 M©. The bolometric light curves were also used to derive 
the evolution in different filters, by artificially setting the object 
to a distance of lOMpc. Finally the same rise time fitting algo¬ 
rithm was applied to the model light curves as was previously 
applied to the SN data. 

In comparison, more sophisticated hydrodynamic al mod¬ 
els allow a proper treatment of the time-dependent radiative 
transfer equations, however they are very time consuming in 
their cal culatio ns, in p articu lar over a larger parameter space. 
iBlinnikov & Bartunovl d!993l) . for example, assume explosion 
energies of E ~ 1-2 xlO 51 erg, ejecta masses of M ~ 1-6 M©, 
progenitor radii of R (l ~ 500-600 R© and derive B-band peak 

magnitudes of-16 mag, which is in good agreement with the 

simple analytic models we use in this paper. 

In Figures [9a] and [9b] we compare our rise time estimates 
with the results from the analytic models. The SNe are depicted 
in the same way as in Figure [8] The shaded regions represent 
the parameter space for specific progenitor radii, respectively. 
For a constant radius the brightest magnitude results from an 
explosion energy of 5 x 10 51 erg and an ejecta mass of 3M©, 
while the faintest magnitude results from an explosion energy of 
0.1 x 10 51 erg and an ejecta mass of 20 M©. 

The models in general predict somewhat longer rise times 
than are observed. It is remarkable, however, that even though 
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Rise time in days (rest frame) Rise time in days (rest frame) 

(a) lArnettl ( j 1980L119821) model (b) iRabinak & Waxmanl J201 lh model 

Fig. 9: Theoretical rise times and end rise absolute magnitudes using the lArnettl (ll980L 1982 .) and iRabinak & Waxman i (1201 it) 
model compared with the SNe in our sample. The SNe are depicted in the same way as in Figure [8] The shaded regions represent 
the parameter space for specific progenitor radii, respectively. The input parameters were varied as follows: 100 Rq < Rq < 900 Rq, 
0.1 x 10 51 erg < E < 5 x 10 51 erg and 3 < M < 20 Mq. The same rise time fitting algorithm was applied to the model optical light 
curves as was previously applied to the SN data. Only those objects that have progenitor detections are labelled explicitly. 


the models used are simplistic in their physical description, the 
results from the models lie reasonably close to the measurements 
from our sample in the magnitude-rise time diagram, and in par¬ 
ticular span similar ranges to the data. The models can therefore 
give us a relatively accurate picture of the trends involved. 

We found that the model parameters M, E, and Rq influ¬ 
ence the peak luminosity and rise time as follows: the luminos¬ 
ity of the optical peak increases with larger Rq and E, and de¬ 
creases with larger M for both the lArnettl (1 19801119821) and the 
IRabinak & Waxmanl (1201 ll) model. For the optical rise time we 
derived approximate analytical expressions that relate it to the 
parameters Rq, E and M (s ee Equ ations IA.111 and IA.12I in the 
Appendix). For the I Arne tt] ( 1980. 119821) m odel we found that 
fopt-rise x w hil e f° r the IRabinak & Waxmanl (120 1 ll) 

model f 0 pt-rise K Rq 55 7’^£j 1 06 M -012 . The dependence of the 
rise time on mass and explosion energy is smaller than the de¬ 
pendence on the radius, as we can see in Equations I A. 1 11 and 
IA.12I In fact our estimate for the rise time shows that to first or¬ 
der the effects of m ass a nd energy on the rise time is negligible 
(in particular for the I Arnet t (; 1 9801 1 1982 ) model). 

With these formulae we now have a quantitative handle on 
the dependency on Rq; it is a square root dependency which fits 
reasonably well with the full solutions (see Fig ures l9al and l9bt. 
Using hydrodynamical simulations ISwartz et al. ( 1991 ) also find 
that larger progenitor radii result in supernovae with longer rise 
times and brighter optical peaks. 

We have seen in the previous section that the difference be¬ 
tween the rise times of SNe II-L and II-P is small, however that 
SNe II-L tend to have somewhat longer rise times and higher 
luminosities than SNe II-P. Additionally, from numerical mod¬ 
els such as iBlinnikov & Bartunovl (119931) we ex pect SN e II-L 
to have higher E /Af values than SNe II-P. In the lArnettl d 198(1 
11 982! ) model a higher E/M gives slightly shorter rise times, so a 


larger Rq is needed to explain the somewhat longer rise times of 
SNe II-L._ 

The IRabinak & Waxmanl 1201 ll ) model, on the other hand 
gives longer rise times for higher E/M values which are of the 
same order (a few days) as the observed differen ce betwe en the 
rise ti mes of SNe II-L and II-P. For the IRabinak & Waxmanl 
(1201 ll) model it is therefore not imminently necessary that the 
progenitors of SNe II-L have larger radii than the ones of SNe 
II-P. 

3.5. Rise times of SNe with direct progenitor detections 

In the following we consider the cases where both direct detec¬ 
tions of a Type II-P/L progenitor and constraints on the rise time 
are available: the two Type II-P SNe 2005cs and 2012A, the Type 
II-L SN 2009kr and the peculiar Type II-P SN 1987A. A sum¬ 
mary of inferred parameters for the progenitors is given in Table 
0 

The parameters can generally be split into two categories: the 
direct observations of magnitudes and colours which were used 
in the literature to estimate the mass, the effective temperature, 
and the luminosity of the progenitor star. We used these in turn 
to calculate the radius of the progenitor star, assuming it radiates 
as a blackbody: L — AnR 2 crT*. The second category consists of 
published models that were tailored to match the observations of 
each SN; this yields estimates of the progenitor mass, explosion 
energy, and the progenitor radius. 

For SN 2005cs the estimates for the progenitor radius vary 
quite significantly, spanning a range from 100 to 700 R© (see Ta¬ 
ble [5] and references therein). It is therefore difficult to compare 
it with the other objects. The two estimates for the progenitor 
radius of t he Type II-L S N 2009kr published bv lElias-Rosa et all 
(1201 Ol) and lFraser et al.l (1201 Ol) are in agreement with each other, 
and it seems that the progenitor of SN 2009kr was larger than 
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that of the Type II-P SN 2012A. iMaund et aD (l2015h however 
argue that the progenitor previously identified as a yellow super¬ 
giant is in fact a compact cluster. Nevertheless, the progenitor 
mass estimates remain within the bounds of previous studies. 
SN 2012A has a relatively short rise time of only a few days. 
The lower limit of the rise time of SN 2009kr is somewhat larger 
than the rise time of SN 2012A, however n ot signi ficantly. 

It is interesting to note that lBlinnikov & Bartunovl ( 1 19931) re¬ 
quire relatively large radii (600-1000 R©) to produce SNe II-L 
with peak fi-band brightnesses between -16 and -17 mag. For 
the relatively bright SN 1979C they even assume a radius of Rq 
- 6000 Rq in their models to match the SN luminosity. 

Even though the results for the progenitor radii are not in 
one-to-one agreement - indeed that was not the intention of 
our stu dy - w ith the values determ i ned va lues using the I Arnett! 
(Il982h or the iRabinak & Waxmanl (1201 lh model, the fact that 
the results are of the same order is encouraging. It is clear that 
a single analytic model cannot fit both cases like a compact pro¬ 
genitor as for SN 1987A, and extended red giants as for some of 
the other SNe. However, Figures l9al and l9bl show that these sim¬ 
ple models succeed in reproducing the parameters of the correct 
order and can help to explain the observed trends as we have also 
seen in Section lA. 3. II 

The results of iGonzalez-Gaitan et al.[ (20151). ba sed on com¬ 
parisons with the analytica l models of Nakar & Sari! d2010h and 
hydrodynamical models of lTominaga et al.l d2009t) . are in agree¬ 
ment with the above, i.e., that the progenitor radii are too small 
when compared with the ex pect ed radii of Type II-P/L SNe. On 
this basis. iGonzalez-Gaitan et all (1201 5l) conclude that temper¬ 


atures of red supergiant stars are underestimated, or that mix¬ 
ing length prescriptions used in stellar evolution models are to 
blame. While this may indeed be the case, we emphasize that 
the very useful, but st ill simp listic framework considered both 
here and in IGonzalez-Gaitan et al . ( 201 5). serve mainly to indi¬ 
cate relevant trends, and that factors such as the explosion energy 
to ejecta mass ratio will strongly influence the actual rise time. 

Interestingly, IGonzalez-Gaitan et al.l d2015l) find that SNe 
with slower declining light curves (i.e. SNe II-P-like) have 
longer rise times than those with faster declining light cruves 
(i.e. SNe II-L-like). This is in direct contrast to our findings. We 
speculate that this difference might arise as a result of inade¬ 
quate constraints on the explosion epoch or photometric sam¬ 
pling of the ris e to p eak for the majority of SNe presented in 
IGonzalez-Gaitan et al.l (l2015i) . 


4. Conclusions 

We presented optical light curves and spectra of the LSQ13cuw, 
a type II-L SN that was discovered unusually early, ~1 d after 
explosion. Deep imaging obtained after the SN had faded away 
revealed a faint extended source, within 0".5 of the reported po¬ 
sition of SDSS J023957.37-083123.8. We propose that this ob¬ 
ject is the host galaxy of LSQ13cuw. 

In order to maximize the amount of information from our 
data, we assembled a sample of 20 II-P/L SNe from the literature 
with well-constrained explosion epochs. We see a tendency for 
SNe that have fainter absolute peak magnitudes (R or r' band) to 
have shorter rise times, and vice versa. 

We performed a simple parameter study using the analyti- 
cal m odels for the light curve cooling p hase from I Arnettl (I1980L 
Il982l) and IRabinak & Waxmaiil (1201 lh . The models in general 
predict somewhat longer rise times than are observed. It is re¬ 
markable, however, that even though the used models are sim¬ 
plistic in their physical description, the results from the models 


lie in a similar region in the magnitude-rise time diagram to the 
measurements from our sample. The models can therefore give 
us a picture of the trends involved. The luminosity of the optical 
peak increases for larger progenitor radii and explosion en ergies , 
and de creases for larger mas s es for both the I Arnettl (I1980L119821) 
and the IRabinak & Waxmanl (1201 lh model. By deriving a simple 
analytical expression to estimate the rise time for both models 
we can show that the dependence of the rise time on mass and 
explosion energy is smaller than the dependence on the progen¬ 
itor radius._ 

In the lArnetld 1980Ul982 ) model a higher E/M gives slightly 
shorter rise times, so a larger Rq is needed to explain the some¬ 
what l onger rise times of SNe II-L. The IRabinak & Waxmanl 
(1201 lh model, on the other hand gives longer rise times for 
higher E/M values which are of the same order (a few days) 
as the observed difference between the rise times of SNe II-L 
and II-P. For the IRabinak & Waxmanl ( 2011 ) model it is there¬ 
fore not imminently necessary that the progenitors of SNe II-L 
have larger radii than the ones of SNe II-P. Independent of the 
model, the early optical light curves show that Type II-L SNe 
have somewhat higher E/M values and/or larger progenitor radii 
than Type II-P SNe. We can rule out that SNe II-L have signifi¬ 
cantly smaller progenitor radii than SNe II-P. They are too bright 
and have too long rise times to fit such a scenario. This is at odds 
with expectations from single star evolution. 
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Table 5: Direct progenitor detections. 


SN 

Type 

SN Prog. 

M in Mq 

from modelling 

E Ro/R© 

[10 50 erg] 

from observations 

T e « L/Lq 

[K] 

R*/R© 


SN 1987A 

II-P pec 

BSG 

initial: 16-22 



15000-18000 

1-041(126 X 10 3 

43 ± 14 

i 

SN 2005cs 

II-P 

RSG 

initial: 9^ 
initial: 10 + 3 
pre SN: 17.3 ± 1 
ejected: 8-13 

56 Ni: 0.0065 
envelope: 9.5 

4.1 ±0.3 

3 

1.6 

600 ± 140 
100 

360 

2800-4000 

1.0-2.8 xlO 4 

200-700 

2 

3 

4 

5 

6 

SN 2009kr 

II-L 

YSG 

YSG 

compact 

cluster 

initial: 18-24 
initial: 15^ 

initial: 16-25 



5300 ± 500 
4300-5300 

1q5.12±0.15 

1q5.10±0.24 

430+ 190 

50o+^o 

-180 

7 

8 

9 

SN 2012A 

II-P 

RSG 

pre SN: 10.5jJ 5 
ejected: 12.5 

4.8 

260 
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Appendix A: Light curve models 

In the following we summarize a few basic concepts that help us 
understand the observed trends in our sample of Type II-P/L SNe 
as discussed in Section^ In this chapter we first have a brief look 
into the different phases of a light curve (Section [Aj} and give 
an overvi ew over two analytic al mo dels for the early light curve 
phase by lArne tt] dl98(l 1 19821) and [Rabinak & Waxmanl (1201 ll) 
(Section [A.2l >. In Section [A.31 finally, we investigate the various 
parameter dependencies in the models. 


Appendix A. 1: Light curve phases 

Radiation hydrodynamic calculations of the expl osio n and e vo 
lution o f Type II SNe have b e en pr e sented bv iGrassberg et al 


1971). lFalk& Arnett 


1983 


FaM (11978b. iKlein & Chevaliet 


1978 l. lffillebrandt & Mullen d 198 lh . and Litvinova & Nadezhin 
1985), as well as more specifically for SN e II-P 


by Youn g! d2004b. llJtrobin d2007h. U trobin & Chugail d2009h . 


iBersten e t alJ (20 1 lk an d iDessart & Hillieri d 201 lh. for SNe II- 
L bv ISwartz et al.1 ( 199lh and Blinnikov &J3artunoy j 1993 ) and 
for the peculiar Type II- P S N 1987Aby BHnnikovetah ( 2000b . 
IDessart & Hillieri (120 1 Oh and lPumo & Zampierild201 lh . Here we 
give a summary view of the physical processes at play and the 
various evolutionary stages seen in these simulations. 

The light curve of a Type II SN can be divided into four 
distinct phases depending on the physical processes at play; 
the shock breakout phase, the cooling phase, the rec omb i natio n 
ph ase, and the nebular phase (e.g. lBlinnikov & Bartunovlll993l : 
lYoungl2004h . The first phase is the shock breakout phase, where 
the internal energy of the shock diffuses out to give an X-ray 
(extended progenitor) or gamma-ray (comp act p rogenitor) burst . 
This phase lasts only for minutes or hours dFalk & Arnettlfl977l) 
and only been observed for two Type II-P SNe (SNLS-0 4D2dc 


and SNLS-06Dlid lGezari et al.1 12008 


and once for a Type I SN (SN 2008D, 


Schawinski et ah 12008b ) 

Soderberg et al.ll2008h. 
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Over the next few weeks, the ejecta expand and cool. As the 
post-shock temperatures are > 10 6 K, hydrogen and helium are 
fully ionized and remain so until the temperature has decreased 
down to ~ 10 4 K, when rec ombination begins. This takes a few 
weeks according to models (iUtrobinll2007l : IBersten et al.]|201 ll) . 
Radiative diffusion is inefficient when compared to adiabatic 
cooling dBersten et al.ll201 H) . 

The ejecta are dynamically evolving during the first few 
days. The adiabatic cooling accelerates the layers from their ini¬ 
tial post-shock velocity to their final coasting velocity, which is 
a factor of v/2 higher (assuming equipartition between internal 
energy and expansion energy in the post-shock flow). The in¬ 
ner regions also experience strong Rayleigh-Taylor mixing that 
breaks spherical symmetry dHerant & Benz! 1199 lb iMiiller et al.l 
11991b . 

The time scale to reach the final dynamic state is approxi¬ 
mately the time scale of adiabatic cooling f ac \ ~ R/v (for an adia- 
batically cooling, radiation dominated gas E ml oc R 1 , where E mt 
is the internal energy) which is 1 day for a radius R - IOOORq 
and a velocity v — 10 4 km s~ ] . Thus, unless the explosion energy 
is very small (giving v <sc 10 4 kms '), it is a good approximation 
to take the ejecta as coasting for t > 1 day. 

Homologous expansion (i: oc R) is achieved on a similar time 
scale, as it is fulfilled when R » initial radius, Rq, the same 
condition for adiabatic acceleration to be complete. 

When the temperature goes below ~ 10 4 K, recombination 
of helium and hydrogen has to be taken into account. By ~ 3 
weeks this begins in the outermost layers of the ejecta and grad¬ 
ually moves into deeper layers as a recombination and cooling 
wave. As the internal energy can easily escape when the gas 
ahead of it recombines and the opacity decreases dramatically, 
energy is now lost by radiation. The recombination wave effec¬ 
tively ploughs through the envelope, releasing internal energy at 
a roughly constant rate giving a plateau in the light curve. 

At some point in the recombination phase, gamma-rays from 
56 Ni located deeper inside the ejecta start to reach the enve¬ 
lope layers and affect the thermal state. In particluar for the later 
parts of the recombination phase, this must be taken into account 
(Eastman et alJ[l994l : lKasen & Wooslevll2009h . 

In order to achieve a linear decline, the recombination 
wave needs to be avoided or weakened. A sufficiently low 
density, for example due to a very low-mass or a very ex¬ 
tended pre-superno va star m ight produce the desired effect 
(Blinnikov & Bartunovll993l) . 

After ~ 100 days the entire ejecta have recombined and be¬ 
come optically thin in the continuum. Internal energy is main¬ 
tained by 56 Co, which heats and ionizes the gas to keep emission 
lines active. 


Appendix A.2.1: Arnett 1980, 1982 

[Arnettji d 19801119821) presented the first analytical models for the 
light curves of supernovae. A80 presented models for explosion- 
powered light curves, whereas A82 extended these to include 
56 Ni (there are however also some clarifications in this pa¬ 
per over the use of the A80 models). This framework assumes 
spherical symmetry, homologous expansion, local thermody¬ 
namic equilibrium, a radiation dominated equation of state, con¬ 
stant opacity, and certain constraints on the density profile. The 
first law of thermodynamics is solved with an Eddington outer 
boundary condition. To obtain an analytic solution, it is postu¬ 
lated that the temperature profile is self-similar. 

The luminosity solution for constant density (labelled “A = 
0” in A8(fl) is 


L(t ) = Lq(M, R<), E, k) exp 


r 0 (M, R 0 , k) \ 2 T h (M, Rq , E) j 

(A.l) 


with (see Eq uations 14, 3 7 and 40 in lArnetdll980l and Equations 
19 and 22 in A rne ttlfl982l) 


Lo(M,R 0 ,E,k ) = 5.2x10' 


,43 


/ Rq \ 

1 E ) 

\ 10 14 cm/ 

\2x 10 51 erg / 


0.4 cm 


V“r) [tt\ er g s-1 ’ ( A - 2 ) 

m 2 g \Mq) 


where 

ti,(M,R 0 ,E ) = 
and 


- 1/2 


R 0 £~ 1/2 M 1/2 s, 


(A.3) 


To(M, R 0 , k) = -^—kMRq's, (A.4) 

pc 

where to is the diffusion time scale at t — 0 (and therefore has no 
dependency on E), r/, is the initial doubling time of the radius, 
P is a parameter that depends on the density profile (A80 shows 
that p « 13.7 can be used for a variety of density profiles), and 
k is the opacity. We have here used E = -p,Mu 2 c (the Arnett 
formulation is r/, = Rq/v sc , where v sc is the scaling velocity). 

For t/Th <sc 1, the luminosity is L(t ) ~ e~' lT ". As ro is of the 
order of 10 4 days, the bolometric light curve is initially almost 
flat. For t/Th » 1 the luminosity declines as e~' R 2T ° Th \ where 
Th is the time required to reach a size of twice the progenitor 
size, which is of the order of 1 day. The effective temperature is 
computed by applying the blackbody formula to L{t) and R(f) = 
Rq + y S c t (see Equations 43 and 44 in A80). 


Appendix A.2: Analytical models for the cooling phase 

Analytical models have the advantage of easy implementa¬ 
tion, allowance for physical understanding, and rapid calcu¬ 
lations over parameter-space. Analytical models for the cool¬ 
ing phase (~ 0,5 - 20 da y s) have been presen t ed by several 


Popov 1993; 


, Nakar & Sarif 

20Hl^_ |Rabinak_& Waxmanl 1201 ll Chatzopoulos et alJ 12012 : 
Nagy et al.l 120141) . In this ar ti cle w e focu s on two particu- 
lar models; those of Arnett! (11980! 1 19821) (A80, A82) and 


iRabinak & Waxman (1 2011 ) (R11). 


Appendix A.2.2: Rabinak & Waxman 2011 

The Rll framework starts with a specification of the progen¬ 
itor density profile, and then computes post-shock density and 
temperature assuming a strong, radiation-dominated shock, and 
a analy tical shock propagation formula from lMatzner & McKee! 
( 1999) . The ejecta are then assumed to evolve adiabatically. 

The luminosity is not found from the diffusion equation as in 
A80, but from determining the photospheric location assuming 


9 A80 also gives solutions for other density profiles, but the differences 
are relatively small and will not be explored further here. 
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a constant k, taking the temperature in this layer to equal the 
effective temperature, and then applying the blackbody formula. 

The resulting luminosity and temperature for a red super¬ 
giant progenitor (density power law index n = 3/2 for efficiently 
convective envelopes) are (see Equations 13 and 14 in R11): 


L(t ) = 8.5 x 10 


,42 


Z70.92 d 
£ 51 ^ 0,13 


$' 27 (M/M 0 ) a! X92 5 


G 016 ergs 


(A.5) 


and 


Tph(t) = 1.6 f 


-0.037 


p0.027nl/4 
^51 ^ 0,13 


(M/Mq) 


0.054^.0.28 5 
*0.34 


G 045 eV, 


(A. 6 ) 


where £51 is the explosion energy in units of 10 51 erg, /?o, 13 
is the progenitor radius is units of 10 I 3 cm, f p is a parame¬ 
ter related to the progenitor density profiles (we use f p = 0 . 1 ; 
ICalzavara & Matznen 120041 find f p = 0.079-0.13 for different 
RSG structures), K 0.34 is the opacity in units of 0.34 cm 2 g"” 1 and 
Is is the time in units of 10 5 s. 


Appendix A.2.3: Physical approximations 

Even though the aforementioned models help us in our under¬ 
standing of the light curve morphology, there are limitations 
that need to be considered. One is the assumption of constant 
opacity. Full ionization and pure electron scattering give k - 
0.4cm 2 g _1 (pure hydrogen) and k - 0.34 cm 2 g ^ 1 (solar hy¬ 
drogen and helium mixture). These are reasonable estimates at 
early times, although free -free absorption may contribute as well 
([Rogers & Iglesiaslll992i) . At later times recombination reduces 
the electron scattering opacity bu t line blocking makes signifi¬ 
cant contributions dKarp et alj|1977l) . 

As mentioned above the ejecta are dynamicaly evolving in 
the cooling phase, which is not considered in these models. The 
outer layers of relevance do, however, only accelerate by a factor 
of V2 = 1.4 so this effect should not be strong. 

Weak or moderate mixing of 56 Ni h as n o influence for Type 
II-P light curves during the first weeks dKasen & Wooslevl2009t) . 
However, if there would be very strong mixing, 56 Ni bullets shot 
into the outer envelope could have an impact on the cooling 
phase. 56 Ni becomes important earlier for more compact pro¬ 
genitors (e.g. for SN 1987A, with Rq - 50 Rn, the 56 Ni-drive n 
second peak began already after a week: [Blinnikov et al.ll200d) . 
This limits the applicability of the models to radii > 100 Rq- 

There is mounting evidence that the SN explosion mech¬ 
anism is asymmetric, a nd that the metal cores obtain highly 
asymmetric shapes (e.g . iKiaer et all 120101 : lLarsson et alJ 120131 : 
iGrefenstette et alJl2Q140 . However, the shocks tends to spherize 
at larger mass coordinate, and the outer layers of interest here 
are lik ely close to spherical symmetry (e.g. (Leonard et al]l2001l 

l2006h . 

A radiation dominated equation of state should be an excel¬ 
lent approximation for these early phases (A82). The R11 mod¬ 
els have more realistic density profiles and photospheric position 
than A80, on the other hand they do not include radiative diffu¬ 
sion which A80 does. 


Appendix A.3: Parameter dependencies in the models 

Here we discuss the influence of model parameters M, E, and Rq 
on light curve properties (bolometric and optical rise times and 
peak luminosities). 



Fig. A.l: Fraction of the of the bolometric light, L opt , that 
emerges in the optical (4000-8000 A). / opt peaks when r p h 0 t ~ 
7000 K (dotted line). 


Appendix A.3.1: Radius 


The dependency of the bolometric luminosity on the progenitor 
radius is L oc R 0 both in the A80 and the R11 model. A larger 
radius leads to higher luminosity because adiabatic losses are 
smaller and more of the explosion energy is retained at any given 
time. 

Both in A80 and R11, larger Rq also leads to brighter optical 
peaks, driven by the dependency of Lboi- A larger Rq additionally 
means it takes longer to rise to optical peak, which counteracts 
this, but not enough to reverse the trend. 

The optical light curve can be thought of as the bolometric 
light curve multiplied by the fraction of the bolometric light that 
emerges in the optical: 


A„pt(0 — £boi(0 x ftp\(t ). 


(A.7) 


The first function we have explicit equations for, the second is 
computed from the photospheric temperature. The optical peak 
will be the peak of / opt (f), shifted to a somewhat earlier time 
by the monotonically decreasing £boi(0 function. Under the as¬ 
sumption of constant opacity, the / opt function depends only on 
the photospheric temperature, which in turn depends on £boi (0 
and R p hot(0- For R » Rq, / op t(0 depends on Rq only through its 
influence on Lboi- 

As exemplified by the Rll models, the bolometric lumi¬ 
nosity follows L(T) ~ tZ 16 and changes on a time scale of 
tl = L/dL/dt - 6.3 t. For our typical rise times of 5-15 days, 
the bolometric luminosity thus changes on a time scale of 32- 
95 days. The temperature instead changes on a faster time scale 
tj = Tph/dTph/dt - 2.2 1 corresponding to 11-33 days for our 
typical rise times. The time scales in the A80 formalism show 
similar relations. This means that the optical light curve peak 
must lie close in time to the peak in / opt . 

We define optical as the 4000-8000 A range. Figure IA.1I 
shows / opt for this range as a function of T. / opt peaks when 
T P hot ~ 7000 K. Indeed, we find that for the A80/A82 model the 
rise to optical peak mostly ends at temperatures between 6400 
and 7200K, and for the Rll model at temperatures between 
7500 and 7900 K for a range of parameters (see Section [T2l for 
the exact parameter space). 

We further assume the peak of L opt is close to the peak of 
/o P t, and R » Rq. Then 


Fphot - 


Abol(0 


4tot (tfphod) 


1/4 


(A. 8 ) 
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where u p h 0 t is the photospheric velocity and <r is the Stefan- 
Boltzmann constant. As Lboi = CRqEM exp (-f 2 /(2roTh)) for 
t » r h ( Arnett 1980 1 and w 2 hot * 10/3EM ', we get for 
t <sc 2r 0 T h 


Tphot(0 


C 1/4 /?q /4 

(^cr) 1/4 ti/ 2 ’ 


(A.9) 


where C - 5.2x 10 11 gem -1 s _1 for k = 0.4 cm 2 g 1 . It thus takes 
longer to reach 7000 K (or any other temperature) for a larger Rq 
because Rq increases the scale of Lb 0 i- The dependency on M 
and E drops out, as Lboi and tr hot both depend in the same way 
on them (oc E/M)\ in this way we can understand why E and M 
are not strong drivers of the optical rise times in the models (see 
Figures[9a]and[9b]in Section[3}. 

We can solve this Equation for the rise time: 


t 


opt-rise 


c i/2 r‘ /2 


and we get 


t 


opt—rise 


9.2 d 


Ro 

100 Rq 


1/2 


/ Tpeak \ 2 

(700017 


(A. 10) 


(A.ll) 


The Rll model gives (solving T phot (t) = 7' pc;lk ): 


t 


opt—rise 


6.8 d 


Ro 


100 Rq 


0.56 


/ Epeak \ 

\7000/fj 


77-0.06 

£ 51 


m r ulz 

10 MqJ 

(A. 12) 


The dependencies are similar to the ones derived for A80; we 
note the weak dependencies on E and M. 

With these formulas we can understand why redder bands 
peak later; they have lower 7' pca k- We now also have a quanti¬ 
tative handle on the dependency on Rq it is a square root de¬ 
pendency which fits reasonably well with the full solutions (see 
Figures [9a] and [9b] in Section [3]). Using hydrodynamical simu¬ 
lations ISwartz et alJ (1199 ll) also find that larger progenitor radii 
result in supernovae with longer rise times and brighter optical 
peaks (/1-band). However, as these authors discuss, significant 
He enrichment of the envelope could have a similar effect to 
smaller envelope masses in that the light curve will peak ear¬ 
lier, be brighter, and still have a linear decline. The reason is that 
helium has a lower opacity than hydrogen and therefore the en¬ 
ergy trapping is less effective. This means that in principle the 
observed brightness for SNe II-L might also be explained by a 
He-rich envelope. In our parameter study we assume a constant 
opacity, k, and therefore cannot account for such effects. 


out more rapidly as higher E leads to higher velocities and lower 
density. On the other hand, higher velocities also lead to stronger 
adiabatic losses, so the final outcome is not obvious. The analytic 
solutions show that the first effects dominate. 

If we ignore the first day or so when the progenitor size af¬ 
fects the time scales, the diffusion time in A80 is 

oc M 3/4 £^ 1/4 . (A. 13) 

Higher mass leads to longer light curves, as it takes longer for 
the energy to diffuse out. Energy has a relatively weak influence, 
but higher energy leads to faster light curves as photon escape 
and adiabatic degrading occur faster. 

Higher mass also leads to dimmer optical peaks (both mod¬ 
els). Higher M scales down Lb 0 i (see above), and this typically 
has a stronger impact than the slower decline rate obtained from 
a longer ToTh- A higher explosion energy, instead, gives brighter 
optical peaks. This is driven by the influence of E on Lboi, as the 
evolution of r p h 0 t only weakly depends on E. 

Regarding rise time, A80 gives somewhat longer rise times 
with higher M and lower E (the weak dependency is absent in 
our approximate formula lA. 11 1) . whereas Rll gives shorter ones. 
One difference between these models is that the photosphere is 
fixed in mass coordinate in A80, but moves inwards in Rll. In 
A80 a lower E/M gives higher roTh which means it takes longer 
to reach r pea k (through the exponential factor in Equation [AT}. 
The dependence of the rise time on mass and explosion energy, 
however, is smaller than the dependence on the radius, as we can 
see in Equation ! A.l 11 In fact our estimate for the rise time shows 
that to first order, the effects of mass and energy on the rise time 
is negligible. 


Appendix B: Additional tables 

Table lB.ll lists the local sequence stars from the SDSS DR9 cata¬ 
logue that were used to calibrate the LT and LCOGT g'r'i '-band 
photometry of LSQ13cuw to the SDSS system (see also Section 

I2T1 . 

Table IB. 21 shows r'/R -band properties of SNe that were not 
included in our II-P/L sample, but are shown in Figures [8] and [9] 
for comparison. 


Appendix A.3.2: Mass and explosion energy 

The dependency of the bolometric luminosity on mass is L oc 
;V/ 1 and L oc M 084 in the A80 and Rll models, respectively. 
A higher mass reduces the peak luminosity as it is more difficult 
for photons to diffuse out of more massive ejecta. ISwartz et al.l 
(i 1991 ) also find that higher mass envelopes result in supernovae 
with dimmer optical peaks (fi-band). 

The dependency of the bolometric luminosity on energy is 
L oc E in the A80 model, and L oc E 092 in the Rll model. A 
higher energy leads to higher peak luminosity for two reasons; 
more internal energy is created in the ejecta, and it can diffuse 
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Table B.l: Optical sequence stars 


Star 

R.A. 

Dec. 

mag 

g' 

error 

mag 

r' 

error 

mag 

i' 

error 

1 

02:39:41.21 

-08:28:16.6 

18.640 

0.008 

17.191 

0.005 

16.404 

0.005 

2 

02:39:43.66 

-08:30:21.2 

18.794 

0.009 

18.505 

0.009 

18.416 

0.011 

3 

02:39:50.81 

-08:28:14.9 

18.112 

0.006 

16.681 

0.005 

15.869 

0.004 

4 

02:39:45.81 

-08:33:40.3 

17.402 

0.005 

17.020 

0.005 

16.879 

0.005 

5 

02:39:55.97 

-08:28:56.6 

17.287 

0.005 

16.312 

0.004 

15.961 

0.004 

6 

02:40:02.10 

-08:29:36.2 

17.280 

0.005 

16.814 

0.005 

16.632 

0.005 

7 

02:40:02.22 

-08:30:06.9 

18.213 

0.007 

17.193 

0.005 

16.825 

0.005 

8 

02:40:02.88 

-08:30:34.8 

18.180 

0.007 

17.160 

0.005 

16.779 

0.005 

9 

02:40:05.63 

-08:28:06.6 

17.945 

0.006 

16.494 

0.004 

15.695 

0.004 

10 

02:39:57.24 

-08:34:40.3 

19.135 

0.011 

18.255 

0.008 

17.893 

0.008 

11 

02:40:08.63 

-08:33:27.7 

16.301 

0.004 

15.425 

0.004 

15.140 

0.004 

12 

02:40:12.25 

-08:32:48.9 

16.952 

0.004 

16.480 

0.004 

16.285 

0.004 

13 

02:40:05.99 

-08:28:47.7 

19.115 

0.010 

17.686 

0.006 

17.053 

0.006 


Coordinates and magnitudes as reported by SDSS. 


Table B.2: r'/R -band properties of comparison SNe 


SN 

Type 

E(B - 
Galactic 

V) 

Host 

End of rise 
absolute magnitude 

Rise time in 
days (rest frame) 


SN 1987A 

II-P pec 

0.066 

- 

-14.65 + 0.29 

10.2+ 1.0 

a,b,c 

SN 2009kr 

II-L 

0.07 

- 

-16.86 + 0.15 

>5.0 

a,b,d,e 

SN 1993J 

lib 

0.187 

- 

-17.87 + 0.44 

22.3 ± 0.9 

a,b,f,g,h,i,j,k,l,m,n 

SN 2008ax 

lib 

0.022 

- 

-17.63 + 0.30 

21.6 + 0.3 

a,b,o,p 

SN 201ldh 

lib 

0.031 

- 

-17.39 + 0.01 

20.0 ± 0.5 

a,b,q 


a) this paper; b) NA S A/IPAC Extragalact ic Database; c) Hamuy et al.|| l988|: df Nakano et al.|| 2009|: el Elias-Rosa et al.ll2010t f) 
Rinero et al. 1993^_g} Corwin et al.l ll99 3l hi CenTurion^etalJll 1 993|_J^jHanzl i et_alJl i 1993|_Q_ Tweedy et al.i l 19931 k) DumontetalJ 


1 9931: DlP ressberg er et al.ll 1 993t mi Lewis et al.ll 1994 : n i lBarbon et al. 11*9951 : oi lPastorello et al.i2008T ~p) iTsvetkov et all 20091 ; q) 

Ergon et al.l2014bt 
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